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Abstract. The compatibility of standard and Bohmian quantum mechanics has 
recently been challenged in the context of two-particle interference, both from a 
theoretical and an experimental point of view. We analyze different setups proposed 
and derive corresponding exact forms for Bohmian equations of motion. The equations 
are then solved numerically, and shown to reproduce standard quantum-mechanical 
results. 
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1. Introduction 

Just about fifty years ago, Bohm proposed a new interpretation of quantum 
mechanics jl,. Although standard quantum mechanics (SQM) is statistical and 
nondeterministic, Bohmian quantum mechanics (BQM) is fully causal. In BQM the 
wave-particle duality is resolved. A quantum object is a particle with well-determined, 
although not accurately known, position and momentum. The quantum object's wave 
characteristics are embodied in a quantum potential that acts on the particle and 
is related to the wavefunction. In a system of n particles with total wavefunction 
\&(ri, . . . , r n , t), the velocity of particle i is taken to be [H El Ej 



In BQM, the probability density of finding particles l,...,n at points r 1 ,...,r n is 
given, as in SQM, by the absolute square of the total wavefunction. The BQM 
probability is subjective, that is, it expresses our ignorance of the particles' true 
positions. The continuity equation satisfied by the probability density implies that 
the latter is consistent with Bohmian trajectories. This means that if Bohmian 
particles are distributed according to the absolute square of the wavefunction at time 
to, then their trajectories will transform their probability distribution precisely as SQM 
predicts. The probability of finding a given particle at a given position is the same 
in both interpretations. Thus, statistical predictions of SQM and BQM should be 
indist inguishable . 

Nevertheless, Ghose jl] as well as Golshani and Akhavan jHl El Ej recently proposed 
experimental setups that allegedly lead to different predictions for SQM and BQM. An 
experiment along these lines has been carried out by Brida et al jH] , and the results were 
interpreted as confirming SQM and contradicting BQM. All proposed experiments are 
based on two-particle interference, and make use of symmetrical arrangements of either 
two or four slits. The disagreement between SQM and BQM is seen in the fact that 
SQM statistically allows some pairs of particles to reach detectors asymmetrically, while 
Bohmian trajectories are claimed to forbid this. These conclusions have been disputed 
by Struyve et al |H] and by one of us (THj. The objection is that references [U El El E] 
make use of unwarranted hypotheses on the initial distribution of the particles' positions. 

In this paper we first examine the double-slit and two- double-slit setups proposed by 
Golshani and Akhavan, showing that in all important aspects the latter in fact reduces 
to the former. Next we obtain the equations of motion exactly and point out general 
properties of the velocities. We then compute Bohmian trajectories numerically. The 
trajectories are used to show explicitly how the agreement between SQM and BQM 
comes about. 
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2. Experimental setups 

2.1. Double- slit setup 

Two-particle interference differs from one-particle interference in that the interference 
pattern does not show up in the individual detection of a particle on a screen but in the 
joint detection of a pair of particles. That is, the interference pattern is a property of 
configuration space. 
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Figure 1. Double-slit setup. Correlated pairs of particles emitted from source S go 
through slits A and B and are detected on a screen. 

The double-slit setup proposed in jU E] is shown in figure HJ A pair of identical 
particles emitted by a source S impinges on a double-slit interferometer. Just behind 
the slits the wavefunction is assumed to be well described by plane waves with zero 
total momentum in the y-direction and identical momenta in the x-direction. Thus 
each particle goes through a different slit. Since both particles are identical, the 
total wavefunction is symmetrical (bosons) or antisymmetrical (fermions) under particle 
permutation. The edges of the slits are assumed smooth enough to avoid treating 
diffraction of wave packets as these go through the apertures |2J p. 177]. Each slit is 
taken to generate a Gaussian wave form in the ^-direction, while leaving the wave form 
unaffected in the x-direction. Explicitly, let particle i emerge form slit A at t — 0. Its 
wavefunction is then written as 

Va(^,2/«* = 0) = (27ra 2 )~ 1/4 exp + i[k xXi + k y ( yi - y)]| . (2) 

We should point out that once we settle on equation (J2J) as the one-particle wavefunction, 
no information on the particle's state prior to t = is relevant that is not contained 
in (j2J. On the other hand, assuming free propagation with time we get for t > 

4>A(xi,yi,t) = (27rof) 1/4 exp {-(yi - Y - hk y t/m) 2 /(Aa o- t )} 

x exp {i [k x Xi + ky(yi - Y - hk y t/(2m)) - hk 2 x t / (2m)\ } , (3) 

where 

f iht \ 
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The two-particle total wavefunction is taken as 

* = N [V>a(zi, Vi, t) ifofa, y 2 , t) ± ip A (x 2} y 2 , t) 4> B (x u y x , t)\ , (5) 
where N is a normalization constant and 

ip B (x,y,t) =i(> A (x,-y,t). (6) 



In (J5J) the plus sign refers to bosons and the minus sign to fermions. Note that the 
wavefunction \l/ is symmetric (or antisymmetric) with respect to a reflection in the 
y = plane. 

2.2. Two- double- slit setup 

A two-double-slit setup was proposed in [H] and is shown in figure A pair of correlated 
particles leaves source S in a state of zero total momentum. The particles therefore either 
go through slits A and B' or through slits B and A'. Accordingly, up to a multiplicative 
factor the total wavefunction was taken in [3] as| 

* = ij}fSxx,yy,t) 4> B >(x 2 ,y 2 ,t) ±i) A (x 2 ,y 2 ,t)i) B '(x 1 ,y 1 ,t) 

+ ^Pa'(x 2 , y 2 , t) ip B (xi, yi, t) ± ip A '(xi, yi, t) ^b(^2, 1/2, t), (7) 

with 

il)Ar{x,y,t) = ij)jL(-x,y,t), (8) 
ifo>(x, y, t) = ip B (-x, y, t) = ip A (-x, -y, t). (9) 
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Figure 2. Two-double-slit setup. Correlated pairs of particles emitted from source S 
go through slits A, B' or B, A' and are detected on screens. 

It is not difficult to check that if the plus sign is picked in (J7J), the global 
wavefunction can be written as 

* = caa[k x (x x - x 2 ))<S>( yi ,y 2 ,t). (10) 

| This wavefunction is symmetric with respect to a reflection in the y — plane, for both bosons and 
fermions. As with (J5J, we could easily make it symmetric for bosons and antisymmetric for fermions, 
by permuting the sign factors of the third and fourth terms. The following argument would then be 
carried out just as easily, except that the final correspondence between equations 1)110 and (|12[) and 
equation JSJ would involve in the latter both the plus and minus signs. 
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If the minus sign is picked instead, the cosine is simply replaced by a sine. Hence the 
xx, x 2 dependence factors out and it is given by a real function. From equation (JTJ), we 
immediately conclude that x\ = = x 2 . That is, the Bohmian particles do not move in 
the x-direction. 

The reason for this is that with the one-particle wavefunctions written as in (JHJ), 
(JHJ), © an d 0, equation ((7j) does not adequately represent the interference situation 
depicted in figure 121 Indeed all one-particle wavefunctions are infinitely extended in 
both the positive and negative x-directions. Clearly, however, equations © and (JEJ) for 
ipA and ip-Q are intended to hold only where x > d, while (JHJ) and Q for i/; A i and ^b' 
hold only where x < —d. The total wavefunction should therefore be written as 

^ = N'[i() A (x 1 ,y 1 ,t) i^ B f(x 2 , y 2 ,t) + ^p A '(x 2 ,y 2 ,t) i; B (x 1 ,y 1 ,t)] (11) 

if Xi > d and x 2 < —d, and as 

^ = ±N'[ip A (x 2 ,y 2 ,t) ip B '(xi,yi,t) + ip A '{xx, Vi,t) iPb(x 2 , y 2 , t)] (12) 

if X\ < —d and x 2 > d. The two expressions of the total wavefunction do not overlap 
in configuration space, whence Bohmian trajectories associated with one are completely 
independent of the other. As it should, (jllj) and (fT2"|) transform into each other (with a 
plus sign for bosons and a minus sign for fermions) under particle permutation. 

It is easy to check that, up to a constant multiplicative factor, equation (JTTJ) can 
be obtained from (jSJ) (with the plus sign) through the substitution x 2 — > — x 2 . From (JTJ) 
we see at once that all components of velocities are the same with both wavefunctions, 
except that x 2 is transformed into — x 2 . We conclude that Bohmian trajectories in the 
two situations are in one-to-one correspondence, with x 2 being reflected in the yz-plsme. 

Exactly the same argument shows that Bohmian trajectories computed with ()12j) 
and (0) are also in one-to-one correspondence, with x\ now being reflected in the yz- 
plane. 

3. SQM and BQM predictions 

In the remainder of this paper we shall discuss the double-slit setup only, since 
results pertaining to the two-double-slit setup can be obtained by straightforward 
transformation. 

Let one pair of particles leave the slits at t — and arrive at detectors at time t 
(both particles have x-momentum equal to hk x ). In SQM, the probability of finding the 
particles on detectors at points yi and y 2 is given by: 

P( yi ,y 2 ,t) = \y(y u y 2 ,t)\ 2 . (13) 

Suppose first that t is such that \a t \ ~ o"o- Then the spreading of the wave packets 
is not very important. Therefore, the probability of finding both particles on the same 
side of the x-axis is very low. The particles will be detected on both sides of the x- 
axis, rather symmetrically if do is much smaller than Y (see figure Q). Moreover, as 
one-particle wavefunctions overlap very little, interference effects will be negligible. 
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Assume now that |er f | 3> <j . This can be obtained either by taking detectors 
further to the right, or by reducing the value of k x . Then wavefunction overlap becomes 
important and interference effects begin to show up. Detection becomes more and more 
asymmetrical, and both members of a pair can even be detected on the same side of the 
x-axis. 

In BQM, particle velocities can be computed from equation (0) and wavefunc- 
tion (|5]l. The calculation somewhat simplifies if k y = 0, which we henceforth assume. 
Substituting (JHJ) and (jBJ) into (jHJ) and discarding multiplicative factors that do not depend 
on 7*1 and r 2 , we find that 

y 2 1 +y 2 2 + 2Y(y 2 -y i y 



* ~ exp {i [k x (xi + x 2 )]} <J exp 
± exp 



y\ + vl 



2Y(y 2 



Aa a t 

yi 



4:a a t 

Making use of we immediately see that 



m 



x 2 . 



(14) 
(15) 

(16) 



The y-components are trickier. Discarding now the x\ and x 2 dependence of the 
wave function and making use of (jlj). we can write 

exp(-J) {exp(-g) ± exp(^)} , 



where 



(17) 
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From ((TJ) we get 



h 



Im 



m 



_df_ 
dyi 



+ 



dg 
dyi 



-exp(-g) ±exp(g) 
exp(-g) ±exp(g) 



Straightforward manipulations and use of trigonometric identities finally yield 

2hYma 2 sin(fata) ± Yh 2 t sinh(2maoa) h 2 ty\ 



2/i 



+ 



{h?t 2 + 4m 2 o"o) [cos (fata) ± cosh(2mo"Qa)] ' h?t 2 
where the upper sign is for bosons, the lower sign for fermions and 
2Ym ( yi - y 2 ) 



4m 2 



o~. 



4 ' 



Similarly, 



a 



2/2 



h 2 f 2 + 4m 2 a 4 



2hYmal sin(fata) ± Yh 2 t sinh(2m<7oa) 



+ 



h 2 ty 2 



(19) 



(20) 



(21) 



(22) 



(23) 



{h 2 t 2 + 4m 2 o"g) [cos(fata) ± cosh(2mcrQa)] ' h 2 t 2 + \m 2 <j§ 

When t is large enough, yi and y 2 are dominated by the last term in (}2~T]) and 
The behaviour of fermions and bosons therefore coincide. 
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It is easy to check that velocity components satisfy the following relations: 

Vi(x 1 ,y 1 - ) X2,V2, ;t) = -y i (x 1 ,-y 1 ;x 2 ,-y 2 ;t) (i = l,2), (24) 
yi(xi,yr, x 2 , y 2 , ; t) = -y 2 {x\, -y 2 ; x 2 , -y x \ t). (25) 
Furthermore, the ^-component of the centre-of-mass coordinate satisfies 

V{t) = 2 = fi¥ + 4mM - (26) 

This is readily integrated as jHj 

{ 2 \ 1/2 
(2T) 

If two particles emerge from the slits at heights symmetrical with respect to the 
x-axis, then y(0) = 0. Equation (|27j) implies that y(£) = for all t, and the particles 
will necessarily be detected symmetrically Likewise if y(0) is very small, specifically, 
much smaller than oo, y(t) should remain small enough so that the two particles will be 
detected almost symmetrically. 

From such observations, Ghose [I] and Golshani and Akhavan 7 j argued that BQM 
could make experimental predictions beyond what SQM allows. One would only have 
to prepare a number of pairs each with \y(0)\ cr , and detect them at time t. We 
have seen that, if \a t \ o~o, SQM predicts highly asymmetrical detection. Yet with y(0) 
small enough, BQM would predict highly symmetrical detection. 

The flaw in the argument was pointed out in (HI E]] , where it was shown that such 
a selection of y(0) values is incompatible with Bohm's assumptions. It is instructive to 
make the argument fully quantitative. In Bohm's theory, the probability distribution 
of particle positions is given by the absolute square of the wavefunction. Making use 
of (J5J), © and (0), we see that at t — the distribution of y-coordinates is given by 

(ky = 0) 

y 2 , 0) = | Af ^vnio 2 )- 1 {F + G ± 2H} , (28) 



where 



( Vl - Y f + (y 2 + Yf 



(y 2 - Yf +_i yi + Yf 
r 2 o 



G = exp<i-^ (go) 



if = exp| ^ 1. (31) 

Straightforward integration shows that the probability is normalized to one if 

iiv|2 =H i±exp (-3)P (32) 

It is easy to check that (yi + y 2 ) = 0. The standard deviation of y(0) then follows 
from an evaluation of {{y% + y 2 f)- All integrations are elementary, and we obtain 

Ay(0) = i v / ((yi + Z/2) 2 ) = ^|. (33) 
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Hence it is not possible, in Bohm's theory, to select initial positions so that |y(0)| <C a . 
Any such selection scheme amounts either to pick a different initial wavefunction, or to 
make assumptions on the distribution of true particle positions different from Bohm's 
original ones. 

4. Trajectories 

To investigate Bohmian trajectories of pairs of particles governed by wavefunction 
we have computed them numerically. A fourth- and fifth-order Runge-Kutta algorithm 
was implemented to solve differential equations (J2~T|) and ()23|) . both in Mathematica and 
through a special purpose program we wrote in C. Results displayed all use m equal to 
the mass of the electron, cr = 1CT 6 m, Y = 5 a and k y = 0. The distance between slits 
and detectors is set at 0.2 m. From (fTBj) we see that the time needed for a particle to 
reach the detector is inversely proportional to k x . 




Figure 3. 25 pairs of trajectories for (a) Hk x /m = 2 x 10 7 m/s and (b) hk x /m — 
2 x 10 6 m/s. 

Figure El shows a number of pair trajectories for two different values of k x , with the 
result that \o~ t \ = 1.16 er in (a) and |er t | = 5.88 a in (b). For each pair, y-coordinates 
were picked randomly according to Gaussian distributions with standard deviation ctq, 
one centred about Y and the other centred about —Y. In (a) the Bohmian trajectories 
are almost horizontal straight lines, which bears out the fact that wave packets have 
spread very little. In (b), on the other hand, Bohmian trajectories spread and a number 
of pairs display considerable interference. Some particles coming out of the upper slit 
end up below the x-axis, and some coming out of the lower slit end up above. 

Figure E] shows a much smaller number of pairs, this time labelled independently. In 
(a) the initial ^/-coordinates of both particles are picked symmetrically, that is, y(Q) = 0. 
Trajectories remain symmetrical, as equation (|2*7J) predicts. In (b), on the other hand, 




Figure 4. Three pairs of trajectories with (a) symmetrical and (b) asymmetrical 
initial conditions. 



the initial y-coordinates of the upper particles are set at y = Y, while the ^-coordinates 
of the lower particles are picked as y = —Y + 1.5 a (solid line), y = —Y (dotted 
line) and y = —Y — 1.5 oo (dotted-dashed line), respectively. The particle leaving from 
y = —Y + 1.5 do clearly ends up above the x-axis, and it gets arbitrarily far from the 
axis if allowed to go on. 

To sum up, the statistical distribution of Bohmian trajectories is fully consistent 
with the predictions of standard quantum mechanics. Moreover, pairs of Bohmian 
particles can be detected on the same side of the setup's symmetry plane, and arbitrarily 
far from the plane. 

In the experiment reported in reference jH] , pairs of photons generated by parametric 
down-conversion were allowed to go through a double-slit setup of the type shown in 
figure ^ A number of pairs were detected on the same side of the symmetry plane. On 
the basis of symmetric photon trajectories obtained in [H], and following the analysis 
of electronlike trajectories carried out in jUHj, this result was interpreted as confirming 
SQM against BQM. 

It is clear that our analysis leads to a different picture. We have shown that nothing 
prevents Bohmian particles from being detected on the same side of the symmetry plane. 
Hence, under the assumption that electronlike trajectories are in that respect relevant 
to photons, the results of [S] should be interpreted as confirming both SQM and BQM. 



5. Conclusion 



In this paper, double-slit and two-double-slit setups for two-particle interference were 
analyzed and shown to be essentially equivalent. Bohmian equations of motion 
for pairs of bosons or fermions were then obtained and numerically integrated. 
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The statistical distribution of Bohmian trajectories turns out to reproduce standard 
quantum-mechanical results. Members of a given pair are not restricted to remain 
on different sides of a symmetry plane. Relevant experimental results are therefore 
consistent with both SQM and BQM. 
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